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SUMMARY 


Calculations, "based on linearized supersonic -flow theory and on the 
additional assumption that the Mach cones yaw with the wing, are made to 
determine the effect of thickness on the lateral force and yawing moment 
of a sideslipping delta wing at supersonic speeds. Generalized equations 
are derived for delta wings with rhombic profile and constant thickness 
ratio. Results indicate that, for wings suitable for supersonic flight, 
the contribution of wing thickneSB to the stability derivatives Cy^ 

and Cnp (rateB of change of lateral-force and yawing-moment coeffi- 
cients with angle of sideslip) is small in comparison with the effects 
that may be expected from a vertical tail. 

INTRODUCTION 

Theoretical investigations of supersonic stability derivatives of 
isolated wings have been confined chiefly to the treatment of mathemati- 
cally thin profiles (flat plates) . Finite thickness will have certain 
effects, some of which can be estimated by consideration of small pertur- 
bations. The purpose of the present paper is to present some calculation 
showing the effect of thickness on the lateral-force and yawing-moment 
derivatives of a sideslipping delta wing at supersonic speeds . 

On the basis of linearized theory, equations are developed for the 
stability derivatives Cy^ and C n „ (derivatives of lateral force and 
yawing moment with respect to sideslip) as functions of thickness ratio, 
aspect ratio, and Mach number for delta wings at zero lift. The airfoil 
sections parallel to the plane of symmetry are rhombuses, that is, 
symmetrical double wedges with maximum thickness at 50 percent chord, and 
are of constant thickness ratio. 

In the analysis the Mach coneB were assumed to yaw with the wing. 
This additional approximation introduces considerable mathematical simpli 
fications into the problem. The assumption of the yawed Mach cones 
raises a question as to the quantitative validity of the resultB. 
Although verification of these quantitative results must await a more 
exact analysis, the qualitative results of the present investigation 
are believed to bring out the salient characteristics and signifi- 
cance of the thickness effects. 
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Generalized curves are presented to permit estimation of the 
derivatives for any specific delta wing. Series of curves giving varia- 
tions of the derivatives with Mach number for various aspect ratios are 
also presented. 


SYMBOLS 


General : 


V 

u 

w 

P 

Ap 

q 

9 

M 



a 

c 


t 

A 

m 

b 

S 

A 


Cartesian coordinates 
flight velocity- 

incremental disturbance velocity parallel to x-axis 
incremental disturbance velocity parallel to z-axis 
density of air 

pressure increment caused by disturbance 
dynamic pressure 

disturbance-velocity potential 
Mach number 


root semichord, measured in plane of wing symmetry 

chord length at spanwise station y, measured in a 
plane parallel to plane of wing symmetry 

maximum thickness of section with chord c at spanwise 
station y 

angle of sweep of leading edge 
slope of wing leading edge (cot a) 
span of wing 
wing area 


aspect ratio 



cot A 
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M- = cot A 



A 

4 



1 


P 

Y 

N 

C Y 


C 


n 


R.P. 


angle of sideslip 
lateral force 
yawing moment 
lateral-force coefficient 

yawing-moment coefficient 



real part 


Stability derivatives: 




3 — > o 



^3 / p ->0 


Subscripts: 

( ) yawing moment contribution (taken about wing apex) 

due to lateral forces multiplied by x moment arms 

( ) yawing moment contribution (taken about wing apex) 

due to longitudinal forces multiplied by y moment 
arms 

3 indicates velocity taken in flight direction, 

except when used for stability-derivative 
notation defined previously 


ANALYSIS 


The analysis is based on the linearized theory of supersonic flow 
and the derivation of equations and application of results are therefore 
subject to the usual limitations and restrictions of the theory. In order 
to minimize the mathematical complexities, all pertinent equations are 
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derived specifically for a delta plan form, with rhombic profile. The 
method is applicable, however, to wings of arbitrary plan form and profile. 


Consider the delta wing shown in figure l(a) . The wing is moving 
in a direction parallel to the axis of symmetry, that is, at zero lift 
and zero sideslip. The incremental disturbance velocity parallel to 
the x-axis, u, is readily obtained by methods outlined in references 1 
to 3* The equation is 


&P -r, p mV 

u = = R.p, 

ox n 


W 1 


■r- 


V Ml - m B 


,2t>2 


cosh -i ?. ,:2 PE 2 j + cosh -i x + 

COBh B|y-mx| COSh B|y + mxj 


2w c 


V\l - 4m 2 B 2 \ 


( codf 1 g : .. S - f a gfr + cosh’ 1 
' B|y - 2m(x-a)| B|y + 2m(x-a)[/ 


( 1 ) 


Wn Vp 

where -i and — for the forward wedge and rearward wedge, respectively, 
are found as follows: 



( 2 ) 


It should be noted that equation (1) is applicable to other Mach line 
configurations as well as for that shown in figure 1. 

Now consider the same wing sideslipping to the right at a small 
angle p. (See fig. 1(b).) The Mach cone is assumed to retain the same 
position relative to the wing for both the straight-flight and the side- 
slipping condition. This assumption follows from consideration of the 
discussion given in reference 4. The expression for u for the side- 
slipping wing is still given by equation (1) except that other values are 

W-i Wp 

required for — and — . For small angles of sideslip these values are 
approximately as follows: 

For the right half -wing. 


- |(l + P tan a) 
y = -|(2+ P tan a) 


(3a) 
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and for the left half-wing, 


W. 


-y = ~(l - 0 tan A) 

W 

j- =* -|(2 - P tanA) 


For the sideslipping condition, equation (l ) becomes 


f*2 = B.P. aZ t 
dx rt c 


(1 + 3 tanA) cosh -l x - mB 2 y 


\ 1 - m 2 B 2 


+ ..(I - P tanA) cosh -l 

Jl - m 2 B 2 B |y + mx| 




B|y - mxT 

2 (2 + 3 tan A ) 
\1 - lua 2 B 2 


(3b) 


cosh" 1 x - a - 2mB " 


X - a - 2mB y 
B |y - 2m(x-a)| 


tgfiAl cosh ' 1 x .-. a + - 2 mB 2 r, 
\ l - Um 2 B 2 B|y + 2m(x-a)| 


(4) 


The disturbance velocity in the flight direction (positive forward) 

may be expressed in terms of components along and perpendicular to the axis 

of wing symmetry, ^2 and -^2. The relation is 
~ " dx dy 


Ur = ^2 - 3 |2 

jX dy 


^3 = ^ " p ^ (5) 

The increment in pressure caused by the disturbance is obtained 

from 

Ap = pVup (6) 

which expresses Bernoulli's law with the approximation of small disturbances. 

From equations (5_) and (6), the pressure coefficient ^ for one 
surface of the sideslipping delta wing is found to be 


Ap = 2/^ . 

q V\dx 



( 7 ) 
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The expression for ^2 is given in equation (4) : the term ^2 ma; 

ox oy 

be obtained in the following manner: 

= A f * dx 0 

8y oy i/_oo dx 


The expression for ^ is therefore 

= R.P. 1 v t 1 - | tanA C osh -1 x + agf 
dy jt c ^ „ r 09 B 11 + 


31 C 2 \ 1 - m 2 B 2 


- + , P _ t a nA co 8 h-l x - mB y + p ^ A ( coeh ~l x _ cosh -l x-a 


2 V 1 - m 2 B 2 


B| mx - y| 


2 - S tan A , -1 x - a + 2mB 


2 \ 1 - 4m B 


+ p) cosh- 1 * - a : ft* 3 

Ji - ij^ E l y ' 2m(l ' al 


- Si cosh- 1 a ± 

ll - B I J + 2mU ' a)| 


The lateral- force and yawing-moment coefficients may then be obtained 
by appropriate integrations over the entire wing area, thus, 
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and 


C - N - = 1 - 

n |pV 2 Sb Sb 


JIf% xfcc^jjf g,d»* 


( 12 ) 


where, for convenience, the moments are taken about the wing apex and are 
considered positive for clockwise rotation. The factor of 2 accounts for 
the equal additive contribution of the upper and lower surfaces. Informa- 
tion pertinent to the limits of integration is shown in figure 2. The 
integral expressions for Cy are presented in appendix A and those for C n 


in appendixes B and C. Appendix B presents the first term on the right- 

hand side of equation (l2) which is designated (C n ) $ appendix C presents 

the second term, which is designated (Cn) v * x 

«/ 


for 

and 


The indicated integrations were performed and the resulting expressions 
the stability derivative Cy^ and for the derivatives (^np) 2 

(pnp) are gi Ten as functions of p in table I. 


For moments taken about an arbitrary point Xo 
may readily be derived as 


an expression for 




(13) 


Values of the stability derivative C n are then obtained by adding the 

P 

previously derived right-hand members of equation (13) • 


RESULTS AND DISCUSSION 


It was pointed out in the introduction that the assumption of yawed 
Mach cones raises a question as to the quantitative validity of the 
results. The small triangular regions of the wing between the assumed 
(yawed) position of the Mach cones and the correct (unyawed) position 
have an area proportional to the sideslip angle 0. It appears that the 
basic pressure for 0 = 0 acting on these small triangular regions may 
make a contribution to the stability derivatives C Y p and C n p that is 

of the same order of magnitude as that produced by the increment in 
pressure due to 0 acting over the entire wing. The quantitative results 
obtained herein must therefore be considered as rough approximations. 

The qualitative results, however, which represent the more important 
consideration of the present investigation, are believed to bring out the 
salient characteristics and significance of the thickness effects. 
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Stability derivative Cy A generalized curve showing the variation 


of A — with (i where n = — - 1 is presented in figure 3 . It 

(t/c ) 2 h 

should be noted that the two discontinuities in slope evident in figure 3 
(also in figs. 4 to 6) occur when the Mach lines are parallel to the 
maximum- thickness line and to the leading edge. For any given values of 
aspect ratio, Mach number, and thickness ratio, the stability deriva- 
tive Cy^ is readily obtained by use of this single curve. Variations 

c Y fi 

of 2 — with Mach number for various aspect ratios are presented in 

(t/c) 2 

figure 4. The contribution of wing thickness to the value of Cy for 

suitable supersonic wings (10 percent thick or less) is found to be small 
in comparison with the value of Cy^ that may be expected from a vertical 

tail. 


Stability derivative C n p.- Generalized curves showing the variations 

0 f A 2 — — ^ and with p are shown in figures 5 (a) and 5(b), 

(t/c) 2 (t/c) 2 

respectively. For any given set of parameters, the stability deriva- 
tive C np is obtained by proper use of figures 4 and 5 and equation (13) • 


with p are shown in figures 5 (a) and 5(b), 


Variations of -4- (moments taken about the --chord point) with Mach 

(t/o) 2 3 

number for various aspect ratios are presented in figure 6. As with Cy^j 
the contribution of wing thickness to C n p for wings of small thickness 

is seen to be small compared with the expected contribution of a vertical 
tail, although for some airplane designs the value of Cn» due to thick- 


tail, although for some airplane designs the value of C n p 
ness may be appreciable when compared with the value of C : 


for the 


complete airplane. 
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CONCLUDING REMARKS 


On the basis of linearized theory, calculations have been made to 
determine the effect of thickness on the lateral force and yawing moment of 
a sideslipping delta wing at supersonic speeds. Results indicate that, for 
wings suitable for supersonic flight, the contribution of wing thickness to 
the values of Cy^ and C n p is small in comparison with the effects that 


may be expected from a vertical tail, although for some airplane designs 
the value of C due to wing thickness may be appreciable when compared 
n 3 


to the value of C- 

n 3 


for the complete airplane . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., October 6, 19k8 
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APPENDIX A 


INTEGRAL EXPRESSIONS FOR C Y 


The integral expressions for Cy for the cases indicated in 
figure 2 are as follows: 


Case I (0 g. n ^ 2 ) : 


na 2 m 2 

4p(t/c) 


C„ = - 


2 Y 


. 1 f 


y+2am 


2am 

-y/i 


A dx dy 


/m 


/ r— r J+2am r2em r^~ 

+ 2 . - ( f 1 ^ f m C dx dy + f [ ^ 

\ 1 - km 2 B 2 \^° ^ By+a yA 1 


C dx dy 


where 


and 


A = cosh -l f - .. Iflftr COBh ’ 1 

B|y - mxj B|y + mxj 


1 1 - a - 2mB 2 y x - a + 2mB g y 

C - cosh B |y _ 2m(x - a) | C ° B B|y + 2m(x - a)| 


Case II nS l): 


na 2 m 2 


4p(t/c) 2 ^ \|i - m 2 B 2 J 0 J y/] 


f 7+2am 

1 2am j 2m 


L 


A dx dy 


where A is the function defined for case I. 
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Case III ((i > 1) ; 


2 2 
ita m 


m 2 B 2 - 1 


2am 

2mB-l rjB 


jt dx dy 


f 2am 
2mB-l 


^ y +2am 2am r y +2am 

2m / 2mB-l I 2m 

it dx dy + / / 

y/m J 0 J yB 


D dx dy 


■p. -1 x - mB^y -1 x + mB y 

D = cos - 1 - — - cos 1 — -l . 

B|y - mxj B|y + mx| 
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APPENDIX B 


INTEGRAL EXPRESSIONS FOR (C n ) x 

Integral expressions for the yawing- moment contribution ( c n) 3 
for the three cases shown in figure 2 are as follows: 


Case I (0 < p. < |): 
ita3m3 


P(t/c)‘ 


1 r^am r 2 

^ ’ ' MTT^ Jo 7 ,a. 


y+2am 

2m 


Ax dx dy 


am 


y+2am 


1-mB / 2m 


y+2am 
r2am f 2m 


\|Tl^ 2V - Jo 


Cx dx dy + 


By+a 


/ am J y/m 
1-mB 


Cx 


where 


A = cosh 


-1 * ~ *3 - 2-, - cosh 
B|y - mx| 


-1 x + mB y 
B|y + mx 


and 


o _ p 

, -1 x - a - 2mEry , _i x - a + 2mB y 

C = cosh x — — ~r - cosh - 1 - — — 7 “T7 

B y - 2m(x - a) B y + 2m(x - a) 


Case II (J < p £ l): 


ita3m,3 ^ ^ = - i /" f ^ Ax dx dy 

p(t/c) 2 ^ n)x ^ 


m c B' 


y+2am 
2m 

y/m 


where A is the function defined for case I. 
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Case III (n > l) : 


\ 1 

___ ^c n ; x = - 


P(t/c) 


2am 

2mB-l f yB 


rtx dx dy 


y+2am 2am y+2am \ 

n 2am 2m 2mB-l ( 2m ^ 

/ «x dx dy + Dx dx dy 

' 2am Jy/m -*0 J yB J 

2mB-l ' 


_ _T x - mB 2 y .ii + mB 2 y 

D = cos — : r - cos — : r 

B y - mx| B|y + mxl 
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APPENDIX C 


INTEGRAL EXPRESSIONS FOR (C n ) 


Integral expressions for the yawing-moment contribution 
for the cases shown in figure 2 are as follows: 


p(t/c) 2 ^ n 'y 


(My - 


1 - m 2 B 2 


y+2am y 

p 2am C 2m P2am F 2a ’ 

/ Ay dx dy - / / Aydxdy 

0 J y/m ^0 J y +2am j 

2m / 


ll - ItmS 2 


(r 2am j- 2a 

/ I Cy dx dy 

J 0 J y+2am 

2m 


y+2am 

j' 2am f 2m 

am J y/m 
1-mB 


Cy dx dy - 


am y+2am 
1-mB T 2m 

) J By+a 


Cy dx dy 


A = cosh' 


-1 x - mB y ,_i x + mB y 

- L — 82 cosh x - £ — 

B|y - mx| B y + mx 


C = cosh 


-1 x - a - 2mB 2 y x - a + 2mB^y 

'■pi > r> Vi -*• — ■ ■■ — 


B|y - 2m(x - a)| 


B|y + 2m(x - a)| 
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Case II 


ita^m 2 

P(t/c) 2 


where A 



is the function defined for case I and 


, 1 - a - 2mB 2 y 

E = cos x -7 

By- 2m(x - a) 


cos“l 


x - a + 2mB 2 y 
B|y + 2m(x - a)| 
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Case III (n > l) : 


rta3m. g 

p(t/c) £ 


( Cn V " vrr-7 


2am 

2mB-l ry B 


jty dx dy 


y+2am 

Z 2am f 2m /- 2a/B r yB r2am r2a 

/ ity dx dy - I / ity dx dy - I / ity dx dy 

2am Jy/m f 2am J y+2am u2a/BJy+2am 

2mB-l 2mB-l 2m 2m 


2am y+2am 2am \ 

r 2mB-l f 2m f 2 mB-1 r2a r 2&/B r2a ^ 

+ / / Dy dx dy - / / Dy dx dy - / / Ity dx dy 

v 0 J yB '-'0 / y+2am ,/ 2am o yB I 

/ 


2 / r a / B r r r 

+ [/ / jtydxdy+l I Eydxdy 

\4m 2 B 2 -l \ Jyfam Jo JyB+a 


a/B r 2a 


2am r 2a 


r cam r 
Ja/B j 


y+2am 

2m 


ity dx dy 


where E is the function defined for case II and 


p P 

^ _n x - mB y _i x + mB y 

D = COB 1 — ; f - COS - 1 - — : * 

B|y - mx| B| y + mx| 
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Information pertinent to integration limits for equations 
presented in appendixes. 




/j. = — 


Figure 3 • - Generalized curve presenting the variation of lateral -force 
derivative with a parameter relating aspect ratio and Mach number. 








